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ABSTRACT 


Experimental data for the charging of an air receiver is 
presented and interpreted in detail. The data indicates a 
substantial departure from the adiabatic behavior. The ex- 
perimental results are used to evaluate existing closed form 
expressions for the thermodynamic state of a gas ina re- 
ceiver. A method for experimentally determining the convec- 
tive heat transfer coefficient is developed, evaluated and 
used in conjunction with these expressions. 

The experimental work was performed from March 1969 
through May 1969 at the Naval Postgraduate School, Monterey, 


California. 
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Pee tN PRODUCTION 


Systems using the charging of a receiver are often em- 
ployed today. A common practice when analyzing a receiver 
being charged is to assume that the process is adiabatic. 
This simplifies the governing equations for the system and 
leads to a fairly straightforward expression for the state of 
the gas in the vessel. During the course of the charging 
process however, a substantial temperature difference between 
the walls of the receiver and the gas in the receiver may 
develop. Often the tank's thermal capacitance is of suffi- 
cient magnitude to permit the extraction of large amounts of 
energy from the gas with only a small change in the tempera- 
ture at the outside surface of the tank walls. This obscures 
the apparent effect of heat transfer, for the absence Of a 
noticeable temperature change at the outside of the receiver 
walls may lend credence to the adiabatic assumption whereas 
actually large amounts of heat are being transferred from the 
gas to the receiver walls. Errors introduced by such an 
assumption may lead to serious problems. For example, if 
heat transfer is neglected, calculations for determining the 
amount of gas needed to charge a vessel av lead to a pre- 
dicted value lower than that actually required. 

In order to aid in the analysis of this problem Reynolds 
[Ref. 1] developed a theory which includes the effects of 
heat transfer for determining the thermodv'.mic state of gas 


in a receiver during the charging process. He ueveloned four 


1] 


closed form solutions for the various magnitudes of the system 
parameters. Using his criteria to determine the proper closed 
form solution, many charging processes may be analyzed with- 
out excessive difficulty. The accuracy of these solutions de- 
pends to a great extent on the accuracy with which the average 
convective heat transfer coefficients can be estimated for the 
charging process. 

In his experimental investigation of the blowdown process 
[Ref. 2] Reynolds tested his theory using an h (average con- 
vective heat transfer coefficient for the entire blowdown 
ee determined by taking the value of his heat transfer 
parameter NTU (number of thermal units, see Appendix A) which 
when used in his theoretical equation of state gave the best 
fit to the experimental data. These h values were on the 
order of those predicted on the basis of an assumption of 
steady state turbulent free convection inside the receiver, 
| a a 
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Therefore Reynolds suggested the use of an h bagedmon abies 
assumption in his solution to the charging process. 

The problem still remains, however, that unless h is de- 
termined by an independent means it is difficult to come to 
any quantitative conclusions as to the accuracy of Reynolds' 
theory. 

The purpose of this study has been to derive a method for 
experimentally determining the average convective heat trans- 


fer coefficient between the gas and the receiver, and to 
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apply these results to three of Reynolds' closed form solu- 
tions in order to evaluate their accuracy. The three closed 
form solutions evaluated are for conditions where the heat 
transfer from the receiver walls to the surrounding medium 
can be neglected. These three cases are the ones most com- 
monly encountered. The fourth of these closed form solutions 
deals with the case where the walls of the receiver are very 
thin and the heat transfer from the receiver walls to the am- 
bient medium must be taken into account. This solution was 
not evaluated due to the difficulties in experimentally deter- 
mining the average convective heat transfer coefficient be- 
tween the walls of the receiver and the surrounding medium. 
Also presented in this paper are representative values for 
the various heat transfer parameters determined in this 
undertaking. 

In the following sections of this thesis the theoretical 
and analytical methods and results are discussed. Following 
a description of the charging system, the theoretical model 
due to Reynolds is reviewed and the methods for the experi- 
mental evaluation of the heat transfer coefficient are devel- 
oped. The experimental results are then discussed in the 
light of their relevance to Reynolds' model and their useful- 
ness in engineering problems. A final section is included in 
order to summarize this work and indicate some areas for 


future investigation. 


1S) 


TT. SORSECTIVES 


Experiments on the charging of a gas receiver were con- 
ducted at varying flow rates and heat transfer environments, 
in order to accomplish the following objectives: 

(1) Evaluate the closed form solutions developed by 

Reynolds to approximate the thermodynamic state 
of a gas during charging. 

(2) Develop, use, and evaluate an expression for ex- 

perimentally determining the average instantaneous 
heat transfer coefficient in a gas receiver being 


charged. 


III. DESCRIPTION OF EXPERIMENTAL APPARATUS 


A. GENERAL DESCRIPTION 

The test apparatus (Figures 1 through 4) consisted of an 
aluminum cylindrical gas receiver with a volume of 1.03 cubic 
feet. The complete physical dimensions of the test apparatus 
are listed in Table 1. The vessel was designed according to 
the ASME Boiler and Pressure Vessel Code for a working pres- 
sure of 200 psig. Its top was removable to allow access to 
the interior of the tank, using an O-Ring flange arrangement 
to insure airtightness when the top was in place. The re- 
ceiver was fitted with a flow metering device mounted in the 
center of the tank top. The other fittings consisted of a 
mounting for a thermocouple probe, a pressure transducer 


mounting, and a bleed valve arrangement. 
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This test receiver was mounted so as to be immersed in a 
larger tank which was either open to the atmosphere or filled 
with ice water depending on the particular run. Air for 
charging the test tank was supplied from two 117 cubic foot 
air vessels charged by an air-cooled compressor. The supply 
pressure to the flow metering device was varied from 180 psid 


er 10 psig. 


B. FLOW METERING 

The mass flow rate for each run was determined by maintain- 
ing a critical flow through one of two nozzles or through an 
orifice. These small diameter (D) devices were machined from 
Stainless steel in accordance with a paper by Grace and Lapple 
[Ref. 3] and the discharge coefficients (C,) given in their 
paper were used in these calculations. The diameter and dis- 
charge coefficients are listed in Table 1. The orifice or 
nozzle was mounted on the tank top ina flange type arrange- 
ment that was sealed by O-Rings and connected to a 1 inch line 
that lead to the supply tanks. The stagnation pressure (Py) 
was read by a local pressure tap and gage arrangement mounted 
just upstream of the flow metering device. The stagnation 
temperature (Ty) was likewise found by using a thermocouple 
mounted in the flow just upstream of the metering device. It 
should be noted here that for velocities of the magnitude oc- 
curring in these tests, it was safe to assume that the differ- 
ence between stagnation and local conditions upstream of the 
orifice was negligible. Thus, assuming adiabatic flow throudh 


the orifice or nozzle, the temperature and pressure values of 


5 


Po and To indicated by the apparatus just described can be 


used in the critical flow equation, 


D 
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The air to the flow metering device was controlled by a quick 


acting Jamesbury ball valve. 


C. TEMPERATURE MEASUREMENT 

The test vessel was fitted with a thermocouple probe con- 
Sisting of four 40-gage copper constantan thermocouples. The 
thermocouple housings were designed so as to measure as 
Closely as possible the local temperature at four equally 
Spaced heights in the receiver (Fig. 1). The thermocouples 
were arranged so that the four could be read in series, and, 
at the same time, any one of the four could be read independ- 
ently. The series signal was recorded on channel one of a 
two channel Moseley strip chart recorder, while the signal of 
the thermocouple being read singly was recorded on a contin- 
uous Brown recorder. According to the manufacturer's speci- 
fications thermocouples of this type are accurate to within 
+1.5°R for the temperature ranges found in these tests. 
Since the thermocouple probe readings were only used qualita- 
tively, no further calibration was made. As was mentioned 
above, the stagnation temperature of the inlet air was read 
from a single thermocouple. The signal of this 40-gage copper 
constantan thermocouple was monitored on a Leeds and Northrup 


recorder. 
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D. PRESSURE MEASUREMENTS 

The pressure in the test receiver was obtained by use of a 
Daystrom Wiancko pressure transducer attached to the receiver 
whose Signal was recorded on channel two of the Moseley strip 
chart recorder. The input voltage (21 volts) to the trans- 
ducer was supplied by a Philco power supply and the transducer 
Output signal was adjusted so as to read 10 psi/in on the 
Moseley recorder. The pressure transducer was calibrated 
Originally using a test gage and later using a dead-weight 
tester. When connected to the recorder the pressure trans- 
ducer registered pressures within + 1 psi of the actual pres- 


sure throughout the range of these tests. 


IV. EXPERIMENTAL PROCEDURE 


The experimental program was divided into the following 

three cases: 

l. In the first group of runs the receiver was charged 
12 times under conditions corresponding to the ad- 
iabatic case in Reynolds' theory. The runs con- 
sisted of three charges at each of four possible 
supply tank pressures; 190, 180, 170 and 160 psig. 
The 3/8 inch knife-edge orifice wa. mounted for 
these tests and the total elapsed time for a run 
was approximately one second. The large tank en- 
Circling the test receiver was open to the atmos- 
phere. Thus the rapid chargings gave NTU values 


which dictated the use of Reynolds' adiabatic 


sig 


model for theoretically determining the thermodvnamic 
state of the gas in the receiver throughout the run 
(See Section V-A-3). 

The second group of runs corresponded to the isother- 
mal case described in Reynolds! theory. The 36 runs 
consisted of 3 runs at each of 12 possible supply 
tank pressures ranging from 180 psig to 10 psig, in 
various intervals. The 1/32 inch nozzle was mounted 
in the flow metering apparatus and the run time var- 
ied from 150 seconds to 75 seconds. Once again the 
large tank surrounding the test tank was open to the 
atmosphere. The resulting very small mass flow rates 
were introduced to correspond to NTU values in the 
region of 7, a lower limit for Reynolds' criteria. 
The last set of runs was made with the test receiver 
completely immersed in ice water contained by the 
large tank. The first 9 runs were made with the 

1/32 inch nozzle in place, three each at suvoply pres- 
sures of 180, 160 and 140 psig. The last 18 runs of 
this case were made with the 1/8 inch nozzle mounted 
in the tank top, three each at supply pressures of 
180, 150, 120, 90, 60 amd 30 psig... The ice bathiwae 
to maintain the receiver walls at a constant temvnera- 
ture and thus present the gas with an isothermal sink 
to correspond to the one described in the analytical 


model. 


is 





The values of the supply pressure listed above represent 
the nominal charging pressure. The actual stagnation pnressure 
just upstream of the flow metering device was recorded for 
each run. 

In making a series of runs, the supply tanks were first 
charged to approximately 190 psig. The system was then iso- 
lated from the compressor and the test receiver was allowed 
to reach equilibrium with its surrounding medium. The quick 
acting valve between the supply tanks and the flow metering 
device was then opened. It was kept open until the pressure 
in the test tank reached a value insufficient to maintain a 
critical pressure ratio across the flow meter. The valve was 
then closed and a run terminated. For the succeeding run the 
supply tank pressure was adjusted to the desired value and 
the test receiver was vented to the atmosphere. The next run 
did not commence until the temperature of the gas in the re- 
ceiver (as indicated by the thermocouple nrobe) reached equi- 
librium with that of the tank walls and thus the surrounding 
medium. Having made one pass through the supnply pressures 
indicated for a particular case, the supply tanks were re- 
charged and the procedure repeated for a second and third 
time. 

Throughout the charging process the pressure and temvera- 
ture history of the gas in the tank was recorded on a two 


channel strip chart recorder. 


dee, 


V. THEORETICAL CONSIDERATIONS 


Aye  LHBORY DEVELOL Bees. REYNOLDS* 


lis) Sine wOOUCELe I 

In an analysis of the charging process one object is 
to obtain the time dependent thermodynamic state of the qas 
in the receiver. A common practice is to express the tempera- 
ture of the gas (the dependent variable) as a function of the 
mass (the independent variable). For a constant volume re- 
ceiver the perfect gas law may then be used to express the 
pressure as a function of the mass. Knowing the mass flow 
rate the thermodynamic state of the gas as a function of time 
can be described. 

With this in mind Reynolds develoned the following 


model for a receiver being charged: 
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The thermal capacitance of the receiver walls and any other 


internal structure is lumped into a single cavacitance 


ithe nomenclature used in this section is taken from Ref- 
erence 1, see Appendix A for listing. 
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represented by Cy and equal to the product of the mass of the 
Capacitance and its specific heat (C = Moo,) The heat trans- 
fer resistance between the receiver walls and the surrounding 
environment is represented by R, and the heat transfer resist- 
ance between the capacitant material and the gas in the re- 
ceiver is represented by R,. 
Z. General Differential Equation €or Charging 

Reynolds assumes that the thermal resistance as well 
as the thermal capacitance are invariant with time. He also 
assumes perfect mixing of the injected gas and the gas in the 
receiver, and that the walls of the receiver are at a uniform 


temperature throughout. Using these assumptions he extends 


his model to include all the significant energy terms as 






follows: 
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Combining an energy balance for the above model and the heat 
transfer rate equations, 
L_=-T T=-T 
es 12 0) 


too “ ~R ego 
58. LE 








the following non-dimensional general differential equation 


for the charging is obtained (see Appendix B): 
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For the case of constant mass flow rate (w* = 1) this equation 


reduces to: 


qa? T* 
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w or + Ga eee =O (1) 
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where Ci Coy C3, and Cy are constants comprised of the sys- 


tem parameters NTU, NTU and C,*; defined as, 
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The parameter NTU is used by Reynolds to represent the conduct- 
ance (reciprocal of resistance) of the receiver as well as to 
give a measure of the rate of the process. The Cy* parameter 
is the ratio of the capacitance of the receiver to the initial 
Capacitance of the gas. 

Reynolds solves equation (1) analytically, however 
the resulting series solution is seen to be impractical for 
many engineering applications. Fortunately in many cases the 


values of the system parameters are such that a simoler closed 
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form solution may be obtained. These solutions can be obtain- 
ed by simplification of the series solution or by returning 
to the general differential equation for charging. In his re- 
port Reynolds nonetheless chooses to develop each closed form 
solution from a simplified form of his original model. A 
brief description of these and their solutions are presented 
in the following four sections. A complete derivation of the 
general differential equation for charging as well as deriva- 
tions of the closed form solutions are presented in Appendix 
Br. 
3. Adiabatic Charging 

As has been mentioned, a common method of ceo se of 
a system being charged is to assume the process is adiabatic. 
This assumption of no heat transfer to or from the gas may be 
found to represent the case of a high mass flow rate charging 
(wy large) process fairly accurately. This may be thought of 
as due to the fact that the gas has not had time to transfer 
a substantial amount of heat to the walls before the charging 
process is terminated. Even if the high mass flow rate proc- 
ess is fairly lengthy, the effect of the heat transfer is 
small compared to the effect of introducing large quantities 
of mass during the charge. Such an assumption may also be 
justified for the case where the gas has been insulated from 
the thermal capacitance of the receiver walls (h = 0). This 
case corresponds to the situation where NTU 1s verv small, 


and the following model applies: 
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The solutions for the dimensionless temperature and pressure 


are, 
kTy-1 
* — xk 
r iste =) ae (.2)) 
and 
p* = KTS (M*-1) ee |: a (3) 


4. Isothermal Charging 

A receiver with a large thermal capacitance (C 
large) and either charged very slowly (Wo small) or having a 
high hA value, will exhibit an isothermal behaviour. The 
normal temperature increase in the gas due to comoression is 
not observed because the energy is extracted from the gas bv 
the thermal capacitance of the receiver before it can become 
Significant. As described above, a system of this type will 


have large values for the system parameters NTU and Cy*, and 


may be represented by the following model: 











The temperature ratio for charging is by definition, 
i = | | 
and therefore, 
p* = M* , 
>. SCharemag at constant Mass Flow with Heat Transfer 
to - Isothermal Sink 
In some systems the thermal capacitance of the wall 
is much larger than the thermal capacitance of the aas (Cy* 
very large). For this case it.can reasonably be assumed that 
the energy extracted from the gas is not great enough to sig- 
nificantly affect the temperature of the receiver walls (TO 
is a constant). Thus the analytical model is modified to 
describe the heat transfer between the gas and an isothermal 


sink as follows: 
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mire selutions for this case are, 


KTF+NTUTS — (kT#-1-NTU+NTUTS) Me EENTO) : 
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ae (4) 


and 


(kD S+NTUT*)M* — (kT*-1-NTU+NTUT*) MeN 
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6. Charging at Constant Mass Flow with Inside Resistance 

Negligible 

A system having a high conductance can be assumed to 
have negligible inside heat transfer resistance (NTU high). 
This implies that the temperature of the capacitance is the 
same as that of the gas and the heat transfer is between the 
Capacitance and the surrounding environment. That is distin- 
guished from the isothermal case in that here the thermal 
Capacitance is finite but not large (very thin walled cylinder 
for example). Since systems of this nature are not too com- 
mon, Reynolds feels that the greatest value of this solution 
is that it supplies information for determining the effect of 
Capacitance on the heat transfer in a charging process. The 


model for this system is: 
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and 
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B. CRITERIA FOR APPLICATION OF CLOSED FORM SOLUTIONS 

In the previous four sections, closed form solutions were 
developed for the state of the gas in a receiver being charged. 
These solutions were developved by simvlifying the original 
model with assumptions as to the magnitude and imvortance of 
various terms in the general solution. For example, in the 
adiabatic charging case the parameter NTU is assumed to be 
very small, thus the mass flow rate is large comwoared to the 
heat transfer coefficient and a solution based on no heat 
transfer is formulated. Similar types of assumptions are 
made for the other three cases. Reynolds then produces quan- 
titative criteria for the use of such assumptions. These 
criteria are based on a maximum deviation of 5% in T* at the 
value of M* = 9. They are determined by comparing the case 
in question with the case which would give the maximum depar- 
ture from this behavior. Thus the adiabatic case. (no heat 
transfer to the walls) is compared to the case where the ca- 
pacitance of the receiver walls is infinite, the isothermal 
sink case (large amounts of heat transfer to the walls). 
When the percentage T* difference between the two solutions 
is plotted against NTU, the value cf the twce solutions aqree 


Vaan ss FOr valtces of NTU Less than .25. In a similar 


fashion criteria are developed for all the closed form solu- 
tions. Some solutions have more than one possible criteria 
for their use, however, only the criteria for the three cases 
tested here which are met by our exverimental system are 
listed below. As mentioned before, the three closed form 
solutions evaluated here are the most useful and the criteria 
for their application are those most commonly found. 
Adiabatic: © =< NIU < 0225 eal C* 


0 


bsothermal Sank: 90225 =a U se Co >a 0 


Tsothermal: NTU > 7 C5 > 40 ‘. 


Cc. DETERMINATION OF CONVEGTIVE HEAT TRANSFER COEFFICIENT 

In order to test Reynolds' theory on the state of a gas 
in a receiver being charged, a technique was devised to obtain 
a value for the convective heat transfer coefficient in the 
receiver in which the tests were run. To this end the follow- 


ing model of the experimental system was constructed: 
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where, 


Ehemmasecw rt low rate into the receiver 


m 

q = the heat transfer to the receiver walls 

d(Mu)/dt = the rate of change of the internal energy 
of the gas in the receiver with respect 
to time 


HS om the stagnation enthalpy of the entering gas 


rH 
il 


the temperature of the gas in the receiver 


P the pressure of the gas in the receiver 


M = the mass of the gas in the receiver . 

Several points should be noted here. First it was assumed 
that the temperature throughout the receiver could be repre- 
sented By a single value T (See ,discussion of imperfect mixing 
section V-E). The same was true for the pressure term P. It 
should also be noted that this development does not include 
the case where there is heat transfer between the surrounding 
medium and the walls of the receiver. 

An energy balance on the receiver yields, 


; 7 _ d(Mu) 
og ae a 


The kinetic energy of the gas in the receiver was neglected 
for it can be shown that for the flow rates in this exneri- 
ment, the velocity of the gas in the receiver was very small 
when compared to the internal energy. Using the notation f 
to represent the derivative of a function f with respect to 


time, the above equation can be written as, 


mH . = Ge= mu + Mu 


Assuming the gas in the receiver is thermallv and calorif- 
ically perfect, 


q = “th(c T-c 7, ) peice (8) 


where T. is the inlet stagnation temperature of the gas. By 
definition of the convective heat transfer coefficient, 


hA(T-T. ) =gqg- “Mc T - m(c T-c T, ) 2 


In the same manner as T represents the average gas temperature, 
T. is used to represent the average temperature of the re- 
ceiver walls. Thus the h defined here represents an average 
instantaneous heat transfer coefficient for all points on the 
inside of the receiver walls. The symbol A represents the 
total inside area of the receiver walls. By use of the per- 


fect gas law it can be seen that, 
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for a constant volume system. Therefore the above equation 


can be written as, 


or, using the perfect gas law again, 


ey ee 
-c P(=)+mhe T. 
4h eee . (9) 


AL (#4) -T 1] 


The next step was to determine an expression for T. se 
was not maintained a constant as in the isothermal sink case. 
In doing so the difficult task of determining exnverimentally 
an accurate average wall temperature was avoided. Since, for 


simplicity, no heat transfer from the outside of the receiver 
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walls to the surrounding medium was assumed, an energy balance 


on the walls yields, 
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where C is the specific heat of the walls and is assumed con- 
stant. Combining this with equation (8), 


q = -m(c T-c T.) - MclT = MCT 


is obtained. Integrating this equation over time from t = 0 


Eo t = © gives, 
(c—9 . a —s ae = Re e 
| ag MCT dt = les oe. m(c_T a Mc T]dt 
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In these tests IT. was a constant, therefore, 
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where subscript 0 denotes the condition at time t = 0. 
Solving for T, gives 
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Again using the perfect gas law to eliminate the term T, the 


final form of the expression for T. 1s obtained, 
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Equations (9) and (10) furnish a means for estimating the 
film coefficient h in terms of the tank pressure and gas mass 
and their rates. The massive experimental simplicity thus 


introduced is the major justification for the acceptance of 


inaccuracies stemming from the definition of the heat transfer 


in terms of spacial averages of gas and wall proverties. 


D. DISCUSSION OF ASSUMPTION OF NO EXTERNAL HEAT TRANSFER 

In the derivation of the expression for h it was assumed 
that there was no heat transfer from the outside of the re- 
ceiver walls to the surrounding medium. This simplifies the 
energy balance, for if there is no heat transfer from the walls 
to the outside medium, all the energy leaving the gas must be 
stored in the receiver walls. This assumption is not only 
valid for this experimental system, but for many pressure ves- 
sels used in charging processes. A vessel capable of with- 
standing a substantial pressure is likely to be constructed 
with materials and dimensions that give it a fairly large 
thermal capacitance when compared to that of the entering 
gas. In the test receiver, for example, the thermal capac- 
itance of the receiver was approximately 8.7 Btu/°R while 
even at its highest mass the gas capacitance was only of the 
order of .08 Btu/°R. Thus the gas temperature drop due to 
heat transfer to the receiver walls would have to be extremely 
large before any significant change in the wall temperature 
was observed. 

This assumption of no external heat transfer may also be 
justified by considering the length of time it would take the 
temperature increase due to heat transfer at the interior sur- 
face to reach the outside walls. One method of analyzing this 
transient conduction heat transfer problem would be to consider 


the receiver walls a semi=—infinite sole with heat etransie: 
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due to convection at the surface. Thus by determining the 
time (t) required for an observable temperature change to 
occur in the semi-infinite solid at a distance (x) from the 
surface equal to the receiver wall thickness of 3/8 inch, some 
idea of the temperature change at the outside surface of the 
test tank will be obtained. According to Carslaw and Jaeger 
[Ref. 4], the solution for this type of problem is given by 


the equation, 
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+ nV 
where 
a = the thermal diffusivity of the material 
h = h/k = the convective heat transfer coefficient/ 
the thermal conductivity 
v = the difference between the wall temperature at 
time t and a distance x from the surface and the 
initial wall temperature 
V = the difference between the gas temperature and the 
initial wall temperature. 
Using the highest value of V recorded throughout our runs 
(100°R) and a corresponding value for h we find that v is 
still only approximately 1.5°R even after 2 minutes have 
elapsed. Thus for our runs the assumption of no external heat 
transfer during a charge is justified. This will be the case 
for all receivers of thermal capacitance sufficient to store 
the thermal energy, transferred from the gas, with negligible 


increase in temperature. 


oS 


E. DISCUSSION OF THRE EFFECTS2OF SEMPRE REEGI sie onic 

In the derivation for determining h, and in the closed 
form solutions derived by Reynolds, it was assumed that the 
effects of imperfect mixing could be neglected. In other 
words, it was assumed that for analytical purposes the gas 
can be considered perfectly mixed so that no temverature or 
pressure gradients exist. As was seen in an unpublished ex- 
periment on the heat transfer in a closed container (with 
Similar dimensions to the test receiver) after gas injection 
[Ref. 5], this is not actually the case. Temperature gra- 
dients do exist and this leads to uncertainty in evaluating 
the heat transfer data correctly. In order to determine the 
effective temperature potential for heat transfer, use must 
be made of some sort of average gas temperature so as to keep 
the complexity of the experimental equipment and the data re- 
duction process reasonable. If this average gas temperature 
is dependent upon the degree of mixing, then the assumption 
that the gas is perfectly mixed may lead to erroneous 
conclusions. 

In his paper [Ref. 1], Reynolds showed through the use of 
a Simple comparison between mixed and ummixed gas systems that 
the pressure is independent of the degree of mixing in an ad- 
jabatic receiver regardless of its magnitude. The magnitude 
of heat transfer in these tests was never very large and, 
Since the average gas temperature used in the derivation for 
determining h is based directly on the perfect gas law, the 
temperature thus calculated is fairly independent of the de- 
Gree Cr iain 
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Reynolds also showed that the temnerature averaged with 
respect to mass in such an adiabatic receiver was independent 
of the degree of mixing. Realizing that this average temner- 
ature would be very difficult to find exverimentally, he went 
on to compare the temperature averaged with respect to mass 
and that averaged with respect to volume and found that the 
discrepancy was so small "that use of a volume average tem- 
perature in experimental investigation is entirely satisfac- 
fory. Therefore, the average temperature measured by the 
four thermocouples connected in series at equal volumes 
vertically in the tank could be used as the temperature refer- 


red to as the gas temperature in the closed form solutions. 


VI. SUMMARY OF EXPERIMENTAL RESULTS 


m™ FORM OF RESULTS 

The data from a particular run consisted of a vressure and 
temperature history for the test receiver during charging, and 
values for the inlet stagnation pressure and temnerature. 
From this data, using the critical flow equation, the mass 
flow rate for each run was calculated and thus the mass of the 
gas could be determined for any instant of time. Using this 
mass and the pressure history of the run, the instantaneous 
h values were determined by equations (9) and (10). The val- 
ues of NTU were then calculated and used to determine which 


class of closed form solution the criteria would designate 


“Re ference ao « 
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for this run. The values of P* = P/P 4 as determined by the 
appropriate closed form solution were then comnared with the 
experimental values. The average percentage deviation be- 
tween these two values and the average value for NTU through- 
out a run were then tabulated. Plots were made comparing the 
P* values as predicted by Reynolds and those observed experi- 
mentally for representative runs. These results are in the 


form of P* versus M* = M/M, because the pressure readings 


0 
were slightly more accurate due to the uncertainties involved 
in determining the temperature of a gas in motion. 

The values of h used in determining NTU have been repre- 
sented by plots of the Nusselt number (NU = AL/k ¢) versus 
eheepmeduce of the Grashor and Prandtl (Gaerne L*p¢g8 -AT/ 

He (Cc u/®) ) numbers for various runs. The subscript f here 
indicates that the term was evaluated at the film temverature, 
which for this case was the mean temperature between that of 


the gas and the wall. The characteristic length L of the 


system was taken to be the receiver height for these 


Calculations. 


B. HIGH MASS FLOW RATE RUNS 

A series of runs was conducted with the test receiver ex- 
posed to the ambient atmosphere at high mass flow rates so as 
to correspond to the adiabatic closed form solution develoved 
by Reynolds. The dimensionless capacitance term 
C6 = C/Mgcys used in the analytical criteria was approx- 
imately 657 for all runs. Based on the criteria NIU - ez. 


for any value of C*, all of these runs were found to correspond 
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to the so called "adiabatic" case. Values of 1.4 for the 
ratio of specific heats of the gas (air in this case) and 1.00 
for the non-dimensional inlet stagnation term Ty, (the inlet 
stagnation temperature of the gas, T, / the original tempvera- 
ture of the gas, Ty) were used in calculating the theoretical 
P* values from equation (3). Figure 5 is a graphical compar- 
ison of the actual and theoretical P* values for representa- 
tive runs of this type. NU versus GRPR plots for a few 
typical runs of this nature are seen in figures ll, 12, 13 
and 14, the points having been plotted at 1/10 second inter- 
vals throughout the runs. Table 2 summarizes the results of 


this series of tests. 


C. LOW MASS FLOW RATE RUNS 

A second group of runs was conducted at low mass flow 
rates so as to correspond to the isothermal closed form solu- 
tion proposed by Reynolds. Once again, the test receiver was 
exposed to the surrounding atmospheric conditions. Reynolds' 
criteria for an isothermal charge solution to hold is NTU > 7 
and C5 > 40. The value of Co for this group of runs was 
again approximately 657 and the NTU values calculated varied 
from 1.8 to 14.6. Thus these runs straddle the cutoff value 
of NTU = 7 and can be used to determine the validity of this 
figure. The solution for an isothermal charge is simply 
P* = M*, for T is a constant and therefore T* = 1. The val- 
ues Of P* predicted by Reynolds and those observed expveri- 
mentally are compared graphically in figure 6 for represent-— 


ative runs of this type. The h values found in a few 
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typical runs are presented by the NU versus GRPR graphs of 
Figure 15, 16, 17 and 18. The points in these figures repre- 
sent 15 second intervals throughout the runs. Table 3 sum- 


marizes the results of this series of tests. 


D. ISOTHERMAL SINK RUNS 

The final group of runs was made, with the test receiver 
immersed in an ice bath, at intermediate mass flow rates so 
as to correspond to Reynolds' isothermal sink model. The 
Criteria for this case is simply Co > 40. However, in order 
to avoid crossing into the adiabatic or isothermal solution 
regions, the value of NTU must be greater than .25 and less 
than 7 respectively. It should be understood that the ice 
bath was not necessary in order to meet Reynolds' criteria: 
with or without it the value of Co (approximately 629 for this 
case) was large enough to dictate the isothermal sink solution. 
The ice bath was employed to ensure that the temperature of 
the receiver walls was a constant and that a significant tem- 
perature potential existed between the receiver walls and the 
gas. Values of 1.4 for the ratio of the specific heats of the 
gas and 1.00 for the dimensionless inlet temperature, Ti, as 
well as for the dimensionless wall temperature, Tar (the re- 
ceiver wall temperature Dey the original temperature of the 
gas Ty) were used in calculating the theoretical P* values 
from equation (5). Figure 7 is a graphical comparison of the 
actual and theoretical P* values for representative runs of 


this type. NU versus GRPR plots for a few typical runs of 


this nature are seen in figure 19, 20 and 21. The points 
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were determined at 15 second intervals for figure 19, and l 
second intervals for the other two figures. Table 4 sum- 


Marizes the results of this series of tests. 


Vit SES CUSSCLON Or RESULTS 


In order to systematically analyze the test results, the 
following section is divided into three separate sections. 
The first section is devoted to an analysis of the equation 
derived for the convective heat transfer coefficient. The 
results obtained when comparing the experimentally observed 
values of P* to those predicted by Reynolds are then discus- 
sed. The last section deals with an analysis of the convec- 
tive heat transfer coefficients obtained with emphasis on the 
identification of general trends that might be investigated 


in a more comprehensive study on this topic alone. 


A. ANALYSIS OF THE EXPRESSION FOR DETERMINING THESCON VEC 
HEAT TRANSFER COEFFICIENT EXPERIMENTALLY 


In an earlier section the expression 


-c B(S) +c T , 
he (9) 


Dp 1 
AL(se)-T 
was derived with the value for Te when it was not considered 
a constant, given by equation (10). 
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As was mentioned before, these expressions were purposely 
developed to obtain a value for h from experimentally deter- 
mined data. Due to the difficulty of measuring the effective 
average temperature of the gas, this term was eliminated in 
the above expressions. The equation for h is thus a function 
of P and M, the other terms being system varameters which were 
generally held constant and known to a high degree of accuracy 
throughout a run. 

First, considering the expression for the wall temperature, 
equation (10), it is seen that even when the temperature of 
the receiver was not held constant the temperature of the walls 
changed very little during a run. This can be seen by noting 
the magnitude of the terms in parentheses as well as the fact 
that the pressure and mass of the gas are never more than one 
order of magnitude larger than their original conditions (Mos 
Py). An investigation of the values calculated for T shows 


0 


that this was indeed the case for these tests. 
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Next the expression for the heat transfer coefficient it- 
eeit wilt be Consiagereda. Since the derivation of this expres= 


sion was centered about the fact that, 


it can be seen that the numerator of expression (9) is an ex- 
pression for the heat transfer q and the denominator is the 
product of the inside area of the receiver walls and the tem- 
perature potential between the walls and the gas. Considering 
first the denominator, it can be seen that when the temvera- 
ture of the gas and that of the walls are the same, this term 
weoezero.. “Of course “in this situation, the heat transfer (the 
numerator of the expression) is also zero and thus as would 
be expected the convective heat transfer coefficient is un- 
defined. Now, examining the case where there is a finite but 
small difference between the wall and gas temperatures, as 
might be the case in a low mass flow rate charging process, 

it is seen that the equation for h becomes very sensitive to 
errors in the temperature of the gas or of the receiver walls. 
It has already been observed that the expression for the tem- 
perature of the walls is insensitive to small errors in the 
values of the system variables, the gas pressure and mass. 
Unfortunately this is not the case for the perfect gas rela- 
tion used to express the gas temperature. Even though an 
error in T may seem negligible when considered on a pnercentaae 
basis, it is not when considering small differences between T 
and T° POusexamolLe, am error of 1°R in the temperature of 


the gas will cut the value for the denominator of h in male 
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if the temperature difference is only 2°R. A temperature error 
of this magnitude is definitely within the exverimental accu- 
racy of these tests, therefore one must exercise caution when 
applying this equation to runs in which the temverature poten- 
tial between the receiver walls and the gas is small. 

The numerator of equation (9) can be analyzed in much the 
Same manner. As was the case with the denominator, it can be 
shown (by examining the normal magnitude of the various terms) 
that except for the case of small temperature differences this 
expression is capable of absorbing small errors in P and M 
without producing a significant error inh. 

Combining these results, it is concluded that equation (9) 
is insensitive to small error in P and M when the temperature 
difference between the gas and the tank walls is at least of 
the order of 10°R or greater. Below this value the calculated 
values of the heat transfer coefficient may be subject to 
substantial error. 

The effect this conclusion has on the confidence in the 
results of these tests is discussed in the following sections. 
The actual value of h is used only in the equation for the 
value of P* in the isothermal sink solution. The runs that 
correspond to this model have substantial temperature poten- 
tials throughout, therefore the values of P* can be calculated 
with confidence. In the other two cases considered, the val- 
ues of P* are not a function of the convective heat transfer 
coefficient, but rather the value of h is used in determining 
the range of NTU in order to determine the particular solution 
to be employed. The values of NTU for these runs are dictated 
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to a greater extent by the mass flow rate than by the convec- 
tive heat transfer coefficient. For this reason, although the 
values of NTU may not be exact, they can certainly be used to 
evaluate the relationship between NTU and the deviation of the 
theoretical P* values from those actually observed in either 
of these cases. 

The situation is less clear cut when it comes to inter- 
preting the convective heat transfer coefficient data. The 
Sensitivity of the expression for h at low values of the tem- 
perature potential places the numerical values for some runs 
in doubt. The low mass flow rate runs in particular are 
greatly affected by this uncertainty. For these runs the tem- 
perature potential was only on the order of 10°R for the ini- 
tial runs and decreased as the charging pressure was progres- 
sively lowered for each successive set of runs. This does not 
mean that this data is not meaningful. There is no reason to 
believe that the trends indicated by such data should not be 
correct; however, it must be realized that the numerical re- 


sults are subject to a large range of experimental uncertainty. 


B. COMPARISON OF REYNOLDS' MODELS WITH THE EXPERIMENTAL 


ReouUnT > 
1. High Mass Flow Rate Runs 


As was mentioned earlier, all runs of this nature were 
meant to correspond to Reynolds' adiabatic solution criteria. 
The value for Cs was 657 in all cases and the calculated val- 
ues of NTU were well below .25 in each run. Reynolds would 


then recommend that the solution to this tvpe of charaing be 
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approximated by the assumption that there was negligible heat 
transfer to the receiver walls. When a comparison of the vai- 
ues of the dimensionless pressure as predicted by this solution 
and those observed experimentally was made the deviation be- 
tween the two values was of the order of 7.5% of the actual 
value (see Table 2), with the theoretical value being higher 
than the actual. 

The discrepancy between these two figures may be at least 
partially accounted for by noting that the assumption that 
there is no heat transfer is not exact. Since the effect of 
heat transfer is to lower the change in internal energy of the 
gas and thus diminish the increase in the gas temperature, a 
lower gas temperature results with a subsequent reduction in 
pressure. With this in mind one expects the theoretical P* 
value to be higher than that actually observed, as is the case 
for this data. 

The parameter NTU is a measure of the thermal conductance 
(h) relative to the rate of gas flow (Wo) « The value of NTU 
should then be related to the exactness of the assumption of 
no heat transfer. A connection between NTU and the deviation 
of the actual and theoretical P* values can be noted ina 
single run. NTU is directly proportional to h for any partic- 
ular rut for Wo is a constant. An increase inh for runs of 
this sort is accompanied by an increase in the deviation (Fig. 
5). Although the quantitative relationshin between h and Wo 
is not clear from the data for different Yuns, the increasaua: 


h and Wo is such that their ratio, and hence NTU, increases 
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with increasing flow. The expected increase in the departure 
of theory from experiment with increasing NTU is observed in 
the comparisons (Fig. 8). 

Reynolds proposed, on theoretical grounds, that for NTU 
less than .25 the deviation between actual values of P* and 
T* and those predicted by the adiabatic solution would be less 
than 5% for values of M* as high as 9. In the experimental 
data, this figure seemed valid at the lowest values of NTU 
(=.08) recorded and appeared slightly optimistic at higher 
values (10% deviation at NTU =.12). In any case, the accu- 
racy of these tests was insufficient to discriminate a 5% 
deviation and the criteria developed by Reynolds is suffi- 
cient for engineering calculations. 

In conclusion, it has been seen that the deviation be- 
tween the actual thermodynamic state of a gas and that pre- 
dicted by the adiabatic solution for charging is a function 
of the actual magnitude of the heat transfer in the system. 
The actual temperature and pressure were always found to be 
lower than that predicted by the adiabatic solution for the 
receiver. The parameter NTU proposed by Reynolds seems to be 
a good measure of whether or not an adiabatic solution should 
be attempted. His criteria seems satisfactory for all but 
the most exacting engineering Poieaiaevens as well. There is 
also evidence to the fact that the values of the convective 
heat transfer coefficient and the mass flow rate are not inde- 


pendent in charging processes of this nature. 
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2. Low Mass Flow Rate Runs 

As was discussed in an earlier section, the runs of this 
type were performed in order to be used in evaluating Reynolds' 
isothermal model. The value for C5 in all cases was 657 and 
the values of NTU ranged from 1.8 to 14.6. For values of NTU 
greater than 7, Reynolds recommends the use of the isothermal 
closed form solution for charging processes. When the values 
of the dimensionless pressure, P*, as predicted by the isother- 
mal solution and those observed experimentally were compared 
the deviation between the two values was somewhere between .5% 
and 4.0% of the actual value (see Table 3). In all cases this 
deviation proved to be negative; that is, the actual pressure 
and temperature were higher than those predicted in an iso- 
thermal solution to the charging process. 

The deviation in the isothermal values of P* and the 
actual values for these runs were small and can be accounted 
for by investigating the assumptions that lead to the iso- 
thermal charging model. In this model the thermal capac- 
itance of the receiver was assumed infinite while the internal 
heat transfer resistance was taken to be zero. Therefore, any 
tendency for the incoming gas to heat up was countered by an 
immediate transfer of heat to the receiver walls. This con- 
dition was approached by charging so slowly that the gas and 
the walls, initially at the same temperature, remained at that 
state. This model can not be exact since the interfaces be- 
tween the gas and the walls have a finite heat transfer resis- 
tance, thus the resulting discrepancies between the actual and 


theoretical values for P*. 
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It seems reasonable to predict in view of the preceding 
discussion that an increase in the heat transfer resistance 
(a decrease in h) would lead to an increased deviation in the 
isothermal and actual values for P*. Although the values for 
h were subject to error in these runs (see section VII-A), a 
definite downward trend throughout a single run was noted 
(Fig. 15, 16, 17 and 18). As predicted, this decrease inh 
was accompanied by an increase in the deviation of the P* 
values (Fig. 6). 

An examination of the effect of NTU on the deviation of 
the isothermal and actual thermodynamic state of the gas was 
then made. It was noted that the values of h for different 
runs remained constant within experimental scatter regardless 
of the value of Wo: Thus, unlike the adiabatic case where 
the effects of changing Wo seemed to be offset by correspond- 
ing changes inh, the value of NTU for this type of charging 
was inversely proportional to the value of Wo: Once again 
NTU proved to be a good measure of the deviation, for as Wo 
decreased the percentage deviation decreased (Fig. 9). 

The criteria proposed by Reynolds suggests that for values 
of NTU larger than 7 an isothermal solution for a charging 
process will give less than 5% deviation from the actual ther- 
modynamic state of the gas. The data from these tests indi- 
cates that for values of NTU as low as 2 the deviation is 
still less than 5%, therefore it seems that Reynolds' criteria 
is slightly conservative. Of course, it must be pointed out 


that the value of Co for these runs was 657 and this is well 
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above the value of Ch > 40 listed along with the NTU > 7 cri- 

teria. Nonetheless, it can safely be asserted that for values 
‘one Ch above 40 and values of NTU of 7 or greater, the isother- 
mal solution to a charging process will give good accuracy. 

In conclusion, it was observed that the actual temperature 
and pressure in a receiver being charged are consistently 
higher than those predicted by an isothermal solution. The 
parameter proposed by Reynolds once again avpeared to be a good 
indicator as to whether or not an assumption of isothermal 
charging can be applied to a particular process. Reynolds' 
criteria also proved to be very satisfactory for engineering 
purposes. Finally, evidence was observed that the value of 
the convective heat transfer coefficient and the mass flow 
rate can be considered independent in a charging process of 


this nature. 


3. Isothermal Sink Runs 
The last group of runs to be evaluated were those made 

with the test receiver immersed in an ice bath. These runs 
were performed at intermediate values of NTU so as to corre- 
spond to Reynolds' isothermal sink closed form solution. The 
only criterion specified for the use of this method is that 
the dimensionless capacitance term be greater than 40. In all 
the runs of this nature the value for Co was 626, therefore 
this criteria was easily satisfied. Of course, Reynolds im- 
plies that the value of NTU dictating an isothermal sink solu- 
tion (to approximate the actual conditions) be be‘%ween .25 and 


7, for if not the regions encountered prescribe either the 
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adiabatic or the isothermal solution. When the values of the 
dimensionless pressure as predicted by the isothermal sink 
solution and those observed experimentally were compared, the 
deviation between the two values was somewhere between 8% and 
14.5% of the experimental value. 

It should be noted here that the isothermal sink solution 
is an intermediate case between the two extremes of an adia- 
batic solution and the isothermal solution. Therefore, the 
thermodynamic state of the gas must be expressed as a func- 
tion of the heat transfer coefficient and the mass flow rate 
in the closed form solution (see equations (4) and (5)). 

As before, physical interpretations were sought for the 
deviation between the isothermal sink solution for the thermo- 
dynamic state of the gas and that actually observed. The 
assumption on which this model was built was that the temper- 
ature of the thermal capacitance remained a constant even 
after absorbing the energy released by the gas in the receiver. 
If the comparatively small amounts of heat transferred to the 
receiver walls were uniformly distributed throughout the very 
large thermal capacitance, this assumption might be quite 
accurate. However, as was seen by an investigation using the 
heat conduction equation for a semi-infinite solid (see sec- 
tion V-D) most of the internal energy (as indicated by the 
temperature distribution) never gets any farther than the 
first few tenths of an inch in a charqing process of this na- 
ture. For this reason the temperature at the wall surface 


does not remain constant but rather increases slightly even 
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though the effective thermal capacitance of the receiver is 
Still extremely large. Bearing this in mind, the actual tem- 
perature potential between the inside surface of the receiver 
walls and the gas is overestimated by the isothermal sink 
solution. Therefore, the amount of heat transferred out of 
the gas and the accompanying temperature and pressure drops 
are slightly exaggerated in this closed form solution (Fig. 
meee 

To begin an evaluation of the parameters that affect the 
deviation between theory and experiment, single runs, in 
which the mass flow rates were constant, were examined. For 
the low mass flow runs there seemed to be only a very gradual 
increase in the value of h as the run proceeds. At higher 
values of the flow rate the increases inh during a run were 
more substantial as were their magnitudes (Figs. 19-21). 
Using both the low mass flow runs (almost constant h) and the 
high flow rate runs (increasing h) it was possible to examine 
the effect of h on the deviation between an actual run and 
the thermodynamic state predicted for this run by the iso- 
thermal sink solution. In analyzing the data (Figs. 7 and 
10) the conclusion was reached that the deviation was appar- 
ently not closely dependent on the value of h. In both the 
low flow rate runs and the high flow rate runs the deviation 
followed the same pattern throughout a run and likewise there 
was no discernible connection between the magnitude of the 
deviation and the values for h. In any case it is probable 


that the effect, if any, of the value of h on the deviation 


of experimental and isothermal sink values is masked by the 
use of NTU = hA/W Cc, aS a parameter in equations (4) and (5). 

For the two extreme solutions (adiabatic and isothermal) 
presented in the previous sections, their accuracy increased 
as the appropriate extreme condition was approached. This 
might lead one to believe that the intermediate ranges of NTU 
for the isothermal sink solution now being considered would 
most closely correlate to the actual thermodynamic state of 
the gas during charging. Since the data does not cover the 
entire range of NTU between .25 and 7, no conclusive statement 
can be made in this regard. Figure 10, however, supports the 
notion that as NTU approaches an intermediate value of say 3.5 
the deviation decreases. One thing can be asserted, and that 
is that, due to the connection between wy and h in this type 
of charging, the values of NTU do not vary greatly. Thus, 
regardless of the connection between NTU and the deviation, 
it seems certain that the deviation between the isothermal 
sink solution and that actually observed is approximately 
constant for a particular system as long as the flow rates 
are of the same order of magnitude. 

In conclusion, it should be noted that the isothermal 
sink closed form solution is a function of the convective 
heat transfer coefficient and the mass flow rate, and thus 
differs from the previous two closed form solutions. The 
deviation between the thermodynamic state of the gas as pre- 
dicted by the analytical solution and that observed experi- 
mentally was between 8% and 15% of the experimental value. 
The closed form solution values were consistently lower than 
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those found in the tests. This discrepancy can be explained 
by the fact that the predicted value of the temperature of 
the thermal capacitance in Reynolds' solution is lower than 
that actually existing. The heat transfer and temperature 
drop are subsequently overestimated in this solution. For 
runs of this nature it was also observed that the values of 

h increase throughout a run and that not only the maaqnitude 
of this increase but the values of h are increased when the 
mass flow rate is increased. There seemed to be little or no 
dependence of the accuracy of the isothermal sink assumption 
on the values for h when Wo is a constant. No conclusive 
statement can be made as to the connection between NTU and 
the deviation of the isothermal sink solution and the actual 
values, although there is evidence to support a prediction 
that this solution is most accurate at intermediate values of 
NTU. Due to the apparent connection between Wy and h the 
values of NTU do not change substantially for runs in this 
region nor do the percentage deviation in the values for the 
thermodynamic state of the gas. Lastly, for this particular 
experiment the discrepancy between the solution as predicted 
by an isothermal sink solution and that observed experimen- 
tally, was slightly higher than that predicted by Reynolds 
even though the value of C* was well above the criteria value 


0 
ene San 


Bie 


ee ee 


Ce SeiSseusstIon OF THE EXPERIMENTALLY OBTAINED HEAT TRANSFER 

COBFPrICcrlENTS 

Having empirically obtained values for the convective heat 
transfer coefficient throughout the various runs, an attempvt 
was made to obtain at least a qualitative explanation of the 
behavior of these values. The charging situation was sub- 
stantially different from those usually described in convection 
theories. The pressure of the system is usually not varied 
in such theories, nor is the direction and magnitude of the 
flow velocity a function of variables other than position. 
Since any attempt to arrive at a theoretical expression for 
the convective heat transfer coefficient in a receiver being 
charged would be completely beyond the scope of this paper, 
no attempt was made. Some insight was gained, however, in 
considering the possible non-dimensional parameters that 
might affect an expression for the convective heat transfer 
in a charging situation. 

From the general equations of motion and energy, an 
average convective heat transfer coefficient evaluated over 
the total surface area of a body (this eliminates any spatial 
dependence) can be expressed as a function of the dimension- 
less parameters as listed below: 

NU = hL/k = £(PR, ER, GR, RE) . 
An investigation of these parameters was then made in order 
to come to a better understanding of the causes for the behav- 
ior of h in the experimental data. First to be considered 
was the Prandtl number. Defined as 
PR = uc /k 


oe 


this dimensionless number can be considered constant and of 
the order of 1 for solutions involving perfect gases. Having 
concluded that the effects of the Prandtl number need not be 
considered, the Eckert number was then considered. This non- 
dimensional parameter expresses the ratio of the inertia terms 
to the enthalpy terms and is defined as 


ER = U*/c AT 


where U is the characteristic velocity of the flow and c AT 
is the change in enthalpy across the boundary layer. In the 
charging processes carried out for this paper the average gas 
velocities were small as compared to the enthalpy, therefore 
the effect of the Eckert number is considered negligible. 

As has been pointed out, a common practice is to take the 
Grashof number as the governing parameter for the heat trans- 
fer processes of the kind under consideration here. Defined 
as, 


GR = gL °BATp*/p? 


this term is the ratio of the buoyancy to the viscous effects. 
The buoyancy of the fluid aids the heat transfer by supporting 
the motion of the fluid at the surface of the receiver walls. 
Noting this, the Grashof number must be considered significant 
in determining the convective heat transfer coefficient ina 
charging process of the sort presented here. Since consider- 
able motion is also present in the tank due to the stirring 
action of the injection, the Reynolds number, defined as 


RE = Uo2/yu 
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must also be considered. In the above expression, U is aqain 
a characteristic velocity of the flow and & a characteristic 
length for the flow. This dimensionless number can be seen 
to be a measure of the relative importance of the inertia and 
viscous effects and often occurs in the analysis of problems 
mn forced” convection. 

Investigating the variables comprising these non-dimen- 
sional parameters leads to some understanding as to their 
magnitude and behavior for charging processes. For small 
temperature changes the density can be considered proportional 
to the pressure. If the temperature changes are significant 
the density is both proportional to the pressure and inversely 
proportional to the temperature. The viscosity is almost 
independent of pressure and increases slightly with tempera- 
ture. For situations of an isothermal nature, it is then ex- 
pected that the variables in the Reynolds number are the 
velocity, the density (which is proportional to the pressure) 
and perhaps, the characteristic length of the fluid flow. If 
the temperature changes throughout a process are large, once 
again the Reynolds number is proportional to the velocity, the 
density (now proportional to the pressure and inversely pro- 
portional to the temperature), the characteristic length, and 
also inversely proportional to the viscosity (which is pronor- 
tional to the temperature). Now applying the same sort of 
reasoning to the Grashof number, the effect of Lt as a variable 
in a moderate temperature change system can again be eliminated 


along with the gravitational constant, gq. The densitv behaves 
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in the manner described in the Reynolds number discussion. 


The coefficient of volumetric expansion, 8, is defined as 


and varies as the inverse of the temperature for a perfect 
gas. The term AT in the Grashof number is the temperature 
change across the boundary layer and L is the characteristic 
length of the system, a constant. 

In preparation for applying what was learned about the 
factors affecting the heat transfer coefficient to the behav- 
ior observed in the experimental runs, the velocity variations 
in the test receiver were examined. When the inlet valve for 
the receiver was first opened a pressure ratio (the receiver 
back pressure/ the stagnation pressure upstream of the flow 
meter) considerably smaller than the critical pressure ratio 
existed across the flow metering device. Sonic flow was 
established at the throat of the inlet to the tank. Therefore, 
fluid entering at a pressure corresponding to the sonic condi- 
tions was forced to adjust to the lower back pressure in the 
tank. Regardless of the form of this adjustment it must have 
been coupled to a substantial increase in the local velocity. 
As the back pressure in the receiver increased, the pressure 
adjustment became less severe and thus the local velocity de- 
creased throughout a charging process. Here the term "local" 
appiies to the conditions in the tank near the point of 
injection. 

For the low mass flow rate runs the temperatv.e potential 


between the gas and the walls, AT, was very small. For this 
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reason it was felt that the Grashof number was not as signif- 
icant as the Reynolds number throughout these runs. Having 
seen that during a charge the velocity near the injector de- 
creased, it was thought that perhaps the characteristic length 
over which the incoming jet acts was also decreased. There- 
fore, although the density increased due to the increase in 
pressure, it seemed reasonable to assume that the Reynolds 
number decreased slightly throughout a run. A drop in the 
Nusselt number, and thus the convective heat transfer coeffi- 
Cient, was predicted. In other words, forced convection gives 
way to free convection as the run slowly proceeds. Figures 
i, 26, ©/ and 18 Seem to support a conclusion of this sort. 
This notion also is supported in a paper by Ulrich et. al. 
[Ref. 6] in which the mode of heat transfer following gas 
injection was observed to proceed from forced to free convec- 
tion. For the high mass flow rate and isothermal sink runs, 
the temperature potential was much greater than that for the 
low mass flow rate runs, thus the Grashof number was expected 
to be more significant in these processes. Since the Grashof 
number is greatly affected by a pressure increase as dictated 
by the p* term, it was felt that the Nusselt number would in- 
crease throughout a run. Once again these assumptions seem 
to be supported by the data, see figures 1l, 12, 13 and 14 as 
well as 19, 20 and 21. 

As to the actual magnitude of the convective heat transfer 
coefficient as indicated by NTU, the data was hard to corre- 
late. As seemed consistent with the conclusions drawn above, 
the values of the Nusselt and Grashof numbers for the low mass 
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flow rate runs were lower than those for the high mass flow 
rate and isothermal sink cases. However, with the exception 
of the set of runs at a very low mass flow rate (perhaps cor- 
responding to the low AT discussion above), the isothermal 
sink values were above those for the high mass flow rate runs. 
One possible explanation for this was the fact that the iso- 
thermal sink runs were recorded in 1 second intervals for run 
times of 7 to 9 seconds. The high mass flow rate runs, on 

the other hand, were only of 1 second duration with the points 
measured at 1/10 second intervals. Therefore it might be that 
for runs of larger duration, the mixing motion caused by the 
impinging jet could spread more thoroughly throughout the re- 
ceiver. These types of forced convection phenomena could not 
be adequately described on the basis of the results of these 
tests. In any case, the values calculated in these tests are 
at most 1 order of magnitude above those predicted by free 
convection theories (Figs. 11-21). Therefore if no better 
estimate of the heat transfer coefficient can be obtained, 

a value based on these theories could be used for processes 
of the nature discussed here. 

The preceding discussion was based mainly on physical 
intuition and is only crudely substantiated by the few experi- 
mental runs performed in this undertaking. As will be 
repeated in the recommendations, a much more extensive and 
controlled study is needed in this area to draw any conclu- 
Sive relationships between the convective heat transfer co- 


efficient and the many variables in a system being charged. 
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The general conclusions to be drawn from this investiga- 


tien are summarized as follows. 


(a) 


(d) 


(e) 


Excellent quantitative agreement between the 
experimental values and those predicted by Reynolds 
[Ref. 1] has been obtained for the thermodynamic 
State of a gas in a receiver being charged. 

The parameter NTU prescribed by Reynolds is a good 
indicator for describing the extent of heat transfer 
in a receiver being charged. 

The technique developed in this paper for experi- 
mentally determining the convective heat transfer 
coefficient may be confidently used in conjunction 
with Reynolds' closed form solutions if the tem- 
perature potential between the receiver walls and 
the charging gas is of an order of 10°R or larger. 
The convective heat transfer coefficient is related 
to the initial mass flow rate of the entering gas 
for all but the extremely slow charging rates 

(Wo <"0i03 lbm/sec.). 

The heat transfer coefficient for charging of the 
nature described in this paper can be estimated 
within one order of magnitude by free convection 
theories. The heat transfer coefficient thus 
obtained will be equal to or less than the actual 


value. The effects of errors in this estimate 
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will be small» provided thatethe: NE Usbses f fine wenn 


small (high injection rates). 


IX. RECOMMENDATIONS 


It is recommended that an extensive and controlled study 
of the effects of the various system parameters on the 
convective heat transfer coefficient during charaging be 
made with emphasis placed on the development of an 
empirical formula for use in this area. 

A possible extension of the technique for experimentally 
determining the convective heat transfer coefficient 
developed in this paper could include the use of analog 
computer methods to directly analyze the pressure and 
temperature signals from probes in the vessel. 

In any further investigation of the type discussed in 
this paper, the following suggestions are made. A more 
exact method of measuring the mass flow rate should be 
developed. Care should also be taken to eliminate any 
possible transient effects in establishing flows used in 
experiments of this nature. It is also recommended that 
hot gases be used in order to obtain substantial temper- 


ature potentials between the receiver walls and the gqas. 
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APPENDIX A 


Nomenclature of Reference l 


English Letter Symbols 


V 


W 


Area, FES 

Thermal Capacitance of receiver shell, Btu/°R 
Specific heat at constant pressure, Btu/(1b°R) 
Specific heat at constant volume, Btu/(1b°R) 
Specific enthalpy, Btu/lb 


Unit heat transfer convective conductance, 
Btu/ (hr £t*°R) 


Thermal conductivity, Btu/(hr ft*°R/ft) 
Mass, lb 

PressUGess Bby £t 7 

Heat transfer rate, Btu/hr 

Absolute temperature, °R 

Total internal energy, Btu 

Specific internal energy, Btu/lb 
Volume, £t° 


Mass flow rate, lb/hr 


Greek Letter Symbols 


8 


Tame, daa 
Vuscositey, Lb/(hr ft) 
Density, lb/ft? 


Denotes a difference 
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Nondimensional Grouping 


* 
CO Co/ (Myc, 
k 2 oe 
* 
M M/M) 
NTU 1/(R, cw 
NTU 1/(R,c Wo) 
* 
P P/P 
* 
ik T/T 9 
* 
ty T/T 9 
* 
aL w/e. 
dts T/T 
w* W/W o 
subscripts 
0 Refers to 
1 Refers to 
S Refers to 
i Refers to 
oo Refers to 


0? 


(hA) . /C wa 


(NA) ../C. Wo 


initial conditions 
inlet state 
Capacitance 

inside receiver 


environmental conditions outside receiver 
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APPENDIX B 


Derivation of General Differential Equation and Closed Form 


Solutions for Charging 


Assumptions and Idealizations 


(1) 


(2) 


(3) 


(4) 


0) 
(6) 


The walls of the receiver are at a uniform temperature 
throughout. 

The heat transfer resistances both inside and outside 
are constant and uniform throughout their respvective 
surface areas. 

Perfect mixing of the injected gas and the gas in the 
receiver. 

The working fluid is a thermally and calorificallv 
perfect gas. 

The specific heat of the capacitance is a constant. 
The kinetic energy of the gas in the receiver may be 


neglected. 


Sememalwitierential @agmation for Charging 







x “sh 
ape 
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ee WAC 
Dd ° ax! 20m 
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An energy balance on the receéiver, using hy as the stagnation 
enthalpy of the inlet gas gives, 


aM _ dau 
OL ee ae meee 


An energy balance on the receiver walls yields, 


The heat transfer rates can be defined in terms of the temner- 
ature potential and the heat transfer resistance as follows: 


Ji, Uy cooly 
Cc 








The change in the total internal energy of the gas in the 
receiver may be written in terms of the mass and the specific 


internal energy as, 


al 92° aa dM 
ici = cic G 


The rate of change of energy in the capacitance may be exvres- 
sed as, 

du aveer.) ar 

ee ee ee 

dé dé Gua 6 


Thus, the energy balance on the gas may be written as, 











T-T 
dM _ ,,du dM @ 
"1 ao - Mae * “ae * Ry 
Or 
T— 
a aM - A eal 4 To 
il; dé dé R. 
and thus, 
T-T 
GMs aT S 
‘epPi-Syt) ag = cvMag * R, 
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Note here that Ty is the stagnation temperature of the inlet 
gas. Therefore, dividing by Cl, gives, 


ar, ee 
dé R. . 
1 


M 





aM 
(kT) -T) a5 


In a similar fashion the energy balance on the canacitance 


can be written as, 











6 lll T-T T -T 
eee, 
aé ie Re F 
Therefore, 
Soe _ ail = ay, i) 
aeé-=6©R.C RC CR .. R . 
sheets. oC oer oo 


The instantaneous mass flow rate, w, may be defined as 
evs GM 
aé 


and we can write, 


as well as, 


cS 
do aM ao “aM : 


The equation for the energy of the gas may then be written as, 





T-T 
_ at 

(kT, -T)w = WM + eae 

i VV 

er 
T-T 
aT SE - 
way + Re pel Sie we Oe CP 


In a similar manner the energy equation of the canacitance 


can be written as 


a ee ee ae 
aM R.C ine Cea ae R 
1 € a Cc ue ol 
Ot 
We S| Le er (2) 
dM (S : R Rec Rise 
Cc alk = nas @- oc 


Expressing eguation (1) in dimensionless form by dividing 


through by To and noting that 


d M d = ee 
“am ~ M mw ~ "ame 
edi) 
M 
0 
we obtain, 
peeps 


* 
wM*S oe + a + (TH-KTH)W = 0. 
Le 





Now dividing by Wo and noting the definition of NTU, we ob- 


tain the dimensionless form, 


* 
w*M#SoS 





+ (T*-T*)NTU + (T*-kTF) w* = 0. (dey) 


In the same way, equation (2) can be non-dimensionalized by 


dividing through by TOW o/Mo so that, 


* * * * 
simul Z acucr, i. ie : 
dM * C war. R Raw Cc Rec we 
Curr re: 00 imo. Cc Eel 


and noting that NTU/CG is egual to My/RWoC, the following 


1s obtained, 
aT * NTU 





wag + TA(NTUSNTU,) - Ae-T* - ae—Te = 0. = (2a) 
0 0 


These two dimensionless equations can be combined to yield a 


single equation giving T* as a function of M* by differentiation 
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of equation (la) with respect to M* and suitable arithmetic 


manipulations represented below: 























* 
weM* STs + (T*-TA)NTU + (T*-KkT#)w* = 0 (La) 
Differentiation with respect to M* yields, 
* 
See ee a att cy 
dmM* dM * dm* dm * ? dm* dm* 
aT dw * dw * 
* ~ ee (¢n = 
‘rearranging, the following is obtained, 
av * a*Ts 
x x ss 
aT * 2 W (250% +m aMez? 4 ar * 
ee NTU dm* 
OL nS Ul ee ae 
. ame ‘2 +M ame ) W kT ame 


NTU ~ NTU 


Substituting this expression for até /dmM* in equation (2a), 





























2mx * * * 
x 2 a = gree aw ae Pucicd 
Sa a a0 ER Cael ae (T*+M* are! 
NL NTU dm* NTU 
dw * 
w*kT* T* NTU 
LamM* NT 00 
- Sar a ax (NTU+NTU,,) a =eT* — AE 1 = Q 
QO 0 0 
Multiplying through by NTU and combining gives, 
aan dw*, aT * 
W*M*ayez + w* (2w*+NTU+M*— >) aoe 
NTU New 
dw* my“ 00 dw * 
+ (w* are - “cx T = Cte = wekT He (3) 


mx 
Cc 
oF ee = 0 
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Now solving (la) for Ts, 


* 
wee + TANTU+T #w*-kT Awe 


dM 
NTU ° 





T* = 
c 
Therefore the term, 


T* 
ee 


(NTU+NTU_)NTU 
oF; co 


in the above expression can be written as, 








NTU+NTU 
aT* ore) TANTU 
* 
W M* ae (ox — cE (NTU+NTU, ) 
kT *w* 


Werle il 
ek: (NTU+NTU_) = —ax- (NTU+NTU,) . 
0 0 





Substituting this in equation (3) we arrive at the general 


differential equation for charging. 


NTU+NTU 








a7*T* dw* 00 aTt* 
kM * * x * 
w*M amM*2 Ny [2w +NTU+M amM* == =Car lage 
_dw* 1 : NTU NTU, . 
+ [w ame + ax (NTU+NTU,.) w + aE: eG 
0 0 
NTU+NTU NTU NTU 
kTs ,dw* Saker k oO mx 
— = ee omer —_ . 
Tiw aM CR kTT _ Ohmme | 0 


Closed Form Solution for Adiabatic Charging 


h, dM 
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Since there are no heat transfer rates to consider in an 
adiabatic solution it is not necessary to specify the mass 
flow rate. An energy balance on the gas in the receiver, 
using stagnation functions, gives 

h,dM = d(mu) = Mdu + udM 


(Gyles 


(A, -u) aM = Mdu 
Then one can write, 


(ea Te=—ceTydM = MexdT 
@ilosé6fYV V 
Dividing by C., gives, 


(kT, -T) dM = Mdt 


In order to get this equation in dimensionless form, we 
divide through by ToMos thus obtaining 


(kT S-T*) dM* = M*dr* 


Now separating variables and integrating this equation from 


fam DCWetial conditions (T = T M=M to a state at a latar 


0) 


0’ 
eime (TI = T, M=M), 

i M_ 

ie "0 aus 

2, OTT ye 

* te 
This gives, 

-log(k,TY-T*) + log(kT}-1) = logM* 


Rearranging, this equation reduces to 


(kT}-1) 
= a * 
7 See 
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Solving for T*, the final result is obtained, 


(kT¥-1) 
a x‘ = 
Ie = kT? —K— 


Thus having derived a relation for T* as a function of M?*, 
the perfect gas relation P* = M*T* is used to obtain, 


p* = kT} (M*-1) + 1 


Closed Form Solution for Isothermal Charging 


Ve 
AFA wi 
T+ T. ° 
_C- 
Since T = constant = Ty = TQ, then 
T# = 7 = 1 
0 


Therefore the isothermal pressure mass relationship can be 
obtained directly from the perfect gas law, and is 


p* = Mx 
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Closed Form Solutions for Charging at Constant Mass Flow with 





Heat Transfer to an Isothermal Sink 


XX 
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An energy balance on the gas uSing stagnation functions gives, 


aM 
1d6 


or Guia). A aM A du, 


h i qo. Wi do do 





By definition of the heat transfer resistance Ray 











nF) 
= c 
Cin are 
al 
and using w = dM/dé, then 
demic du 
Sooo = R, ate Map Tt Woot 
er 
Jeol 
f Cc ‘euk 
S T) Wo. = =a te Mc.a6 + Woe,,t 


di 
Using the following relation, 


at aT dM aT 


de dm do  “oam 
and dividing by Wo and rearranging the exvression becomes, 


T-T 
R;Wo 





OM 
asl + Cc. Man + ct + 


Jak 


Now dividing by Cc. to obtain 


T-T 
dT C 
Mag + tT kT, + EW =~ Q 
re (6) Marg 
Note that all functions refer to stagnation conditions. In 


order to get this equation in dimensionless form it is divided 


by To so that, 
Mas © et * *_T* it 
7 aM + T kTY + (T LS) 
0 ni 0) RY 
Since 
a (=) 
M dT M 0 aT* iL 
a = Se = M*Se and NTU = ——— 
T, MM, ay dM ReWoc,, 
0 


this expression can be written as, 


atT* 
(eed coal Kok kok = 
Mt + (D*=kT*) Gee(T*-T*) NTU = 0 
In order to obtain the closed form solution the variables 
are separated and the expression integrated from the initial 


Condi tionss() = Tos M = My) to a state at a later time (T = 


M = M). 
a M_ 
‘a0 aT* — | ™o ame 
[T* (1+NTU) -kT¥-NTUTS] Se M* 
T . M 
_o _O 
To Mo 
This gives 
a oe = * * 
TLINTU) log | kTS NTUTC+ (1+NTU) T* ] 
Hf —_— 
= TI+NTU) log [-kTF-NTUTA+ (1+NTU) ] = -logM* 


Ta 


@r 


-kT# - NTUT* + (1+NTU)T* = mx” (1+NTU) 


= = * 
7 [ kT, NTUTA+ (1+NTU) ] 


Seiving thas expression for T* as a function of ™* the folow-— 


ing is obtained: 


kTS+NTUTA- (kT 4NTUT*-1-NTU) Me CP FNTU) 


c= 
dies (1+NTU) 
Using the perfect gas law, 


(KD S4+NTUT*) M*— (kT$-1-NTU+NTUT*) M* 


= 
2 1+NTU 


Closed Porm Selutien fer Charging at Constant Mass Plow wath 


Inside Resistance Negligible 





du _ d(Mu) 
dt =e 
OUV 
hom uae 
PE Men 





An energy balance on the gas and the capacitance combined 


yields, 


Se re ae 


h dM d (Mu) ar “ 
a6 dd cdo Teo 


By definition of the heat transfer resistance between the 


Capacitance and the surrounding medium, R,, 








US 


Now since Wy = adM/dé@ the expression 


at _ aT aM _ at 
Bm oct 


can be written. Using this expression and the perfect aas 


relationships in the original energy balance, 


Ave 
dM dM Gull dT 00 
Cpl iae + aecy? + Mevag + Coase + “RO = 


Dividing througihieb, c,,dM/dé or C Wy and noting that AT ./ dé 


= dT/dé, 
S hac 
each co re 
“kT, see hae (iS py aM Ga = - la ane 0 ° 
V vio 0 


In order to express this equation in terms of non-dimensional 


temperature, T*, it is divided by T, to obtain, 


0 


or + cxele 








M* - kT¥ + T*(T*-T*)NTU_ = 0 


OdmM* 1 
The variables are now separated and the expression integrated 
from the initial conditions (T = Tos Mz= My) to a state ata 


later time (i—="T, Me=9M) . 


Th M_ 
Al: M 
ts avs —( 2 ame 
T*(LENTU)-kTFoNTU_T* 7 MFaCe 
S| CBN CEI ae ues A REC 
_0 on’) 
T9 Mo 
The solution to the above gives, 
(M*+C4) 
* ee oe * _- — 
log [T* (1+NTU_) kTy NTOT) ND ee TECR) 


+ log [-kT}-NTU,,TH+ (1+NTU,) | 
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er 


T*(1+NTU,) - kT* - NTU, T* 
une [-kTS-NTU,TSt+ (L+NTU) } 
0 


Solving this expression for T* gives 





+ kT* + NTU T* 
Ali 00 60 





This expression may be combined with the perfect gas law in 
order to obtain an expression for P* as a function of M*, 


L+NTU 
ce 





+ KTTINTULTS 
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Table l 


Physical Dimensions of Experimental Apparatus 


1c 


Test Receiver 

Volume: 1.03 cubic feet 

Internal Area: 5.61 square feet 

Nominal Internal Diameter: 10.19 inches 
Nominal Outside Diameter: 10.75 inches 
Nominal Internal Height: 21.75 inches 
Mass of the Receiver: 38.5 lbms. 
Material: 6061-T6 Aluminum 

Specific Heat: .226 Btu/lbm. °F 

Thermal Conductivity: asl2o Seu. feeecer 
Flow Metering Devices 

3/8 in. knife-edge orifice - discharge coeff. = .620 
1/8 in. nozzle - discharge coeff. = .858 


1/32 in. nozzle - discharge coeff. = .856 


76 


Table 2 


Summary of High Flow Rate Runs 


a: 


ca —s — 
Co = 657 ee TF = 1. © 
Total Run Mass Flow Average Average % P* 
Time Rate NTU Deviation 
(sec) (lbm/sec) 
1 nol 22 bos 
1 325 092 Sodeh 
1. 8825 104 S50 
al - oL0 EOL 7 aw 
1, .309 089 6 os 
1 . 309 085 Gd, 
z . 294 1071, Ld. "OG 
1 {25% 084 5 Se 
1 293 .084 4.92 
1 - 28 162 6.93 
i Bye ~ LO5 7. 3 
1 LO «LES 9.63 


Table 3 


Summary of Low Flow Rate Runs 


Ch = 657 
Run Total Run Mass Flow Average Average % P* 
No. Time Rate NTU Deviationt 
_ (sec) (lbm/sec) 
ut LOS 200s 1 Bn 09 -2.61 
2 LOS OOo: 200 ~3.41 
3 OS 70031 Zao -3.18 
4 LOS .0029 356 9 -2.84 
5 105 0029 Sg Ud =-3.11 
6 FOS 0029 4.78 -2.71 
i IOVS 0028 1.86 =o ne 
8 IOUS: 0028 3.49 -2.93 
9 OS 0028 Ie -3.39 
10 NOS 0024 2a -2.76 
il 1h LOS 0UZs Bh oH, -2.65 
12 105 0025 5.54 -2.22 
diss HOS ~U0Z2Z S234 -2.50 
14 O's 0022 Te7 8 =o 
als) 105 0022 5 6 | -2.09 
16 105 0018 B52 =2202 
17 105 0018 7.80 -1.67 
18 Ubishs 0019 4.43 = 2 ele 
19 105 0015 Seely = 17 0 
20 LOS 0015 8.00 -1.44 
2A. 105 Ciors Cosh -1.27 
22 1205 0012 4.89 -1.54 
23 105 0012 6,02 -1.33 
24 105 0012 S205 orl le 0, 
25 105 0009 4.32 -1.43 
26 EOS VOOs 600i -l1.12 
27 105 .0009 6220 -1.07 
28 105 0007 4.58 -1.49 
Te: TOS 0007 Sd -1.36 
30 105 0007 3.47 -1.45 
Shak 90 0006 4.45 =1.07 
a2 90 0006 ono -0.46 
33 90 0006 iat) -0.52 
34 Fes -0004 14.49 -0.54 
35 WS 0004 14.62 -0.54 
36 1S -0004 259 -0.51 


TR negative deviation indicates that the P* values predicted 


by Reynolds were lower than those actually observed. 
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Table g4 


Summary of Isothermal Sink Runs 


Run Total Run Mass Flow Average Average % P* 
No. Time Rate NTU Deviationt 
“a (sec) (lbm/sec) 
ll IsS50 -AbLO Sil 2.69 -10.39 
2 150 2003 & Dera - 8.84 
B 150 .0031 Zl f -11.08 
4 1650) Ba2 3 A AOY <i . 59 
5 150 Soni 3 2a -10.35 
6 150 2002 & 1 39"7 -11.47 
7 gould BOOZ 5 8.02 = 2 10 
8 150 BOOZ 5 32.30 - 8.58 
9 isye) mOO26 We Sh: - 8.94 
10 9 a0 48:6 ale: = 252 
11 9 m0 4 78 lean -13.40 
Le 9 .0486 ioc — 18.22 
3 9 .0420 hs &S, =1Z67 
14 9 ~0422 lige 8 =I 33 / 
is 9 .0422 iligetdl 5 -14,34 
16 9 20650 1 441 -12.49 
Ine 9 A0) S148 130 3 -13.68 
18 9 F0o49 4 =13.57 
19 9 0274 13 =m 7 8 
20 9 7027 2 1 a3 -12.36 
ZI 9 .0272 142.0 -12.47 
D2 9 -0194 1 35 - 9.92 
23 8 wOur9 2 1.09 =-11.14 
24 8 0190 1 305 -11.14 
Tess y ao a 2 - 7.68 
26 a, sod) 1 me 20, - 7.85 
27 i Peed eS - 7.98 


an negative deviation indicates that the P* values nredicted 
by Reynolds were lower than those actually observed. 
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